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Abstract 



A new procedure of trial variational wave functional is proposed for investi- 
gating the mass renormailzation and the local structure of the ground state 
of a one-dimensional quantum sine-Gordon model with linear spatial modu- 
lation, whose ground state differs from that without modulation. The phase 
diagram obtained in parameters (aA -2 , (3 2 ) plane shows that the vertical part 
of the boundary between soliton lattice phase and incommensurate (IC) phase 
with vanishing gap sticks at (3 2 = Air, the IC phase can only appear for /3 2 
> Air and the IC phase regime is enlarged with increasing spatial modulation 
in the case of definite parameter aA -2 . The transition is of the continuous 
type on the vertical part of the boundary, while it is of the first order on the 
boundary for 1 > An. 

PACS: 05.70.Jk, 68.10.-m,87.22.Bt 
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I. INTRODUCTION 

Recently, many experiments 0~i have detected there exists a soliton lattice phase be- 
tween spatially homogenous commensurate phase with finite gap and incommensurate (IC) 
phase with vanishing gap (free massless boson field) when some quasi-one-dimensional (ID) 
spin-Peierls (SP) materials were placed in an external magnetic field. These experiments 
confirm the 20-year-old prediction of a field-induced transition from the SP state into a spin- 
lattice modulated phase @ and demonstrated that this novel phase has the form of a soliton 
lattice, il The density of solitons can be tuned continuously by adjusting the magnetic field, 
so that their internal structure can be systematically studied. 

In deed, physical realization of SP systems offers the unique possibility to check the 
understanding of this fundamental ID Hamiltonian for arbitrary band filling. Applying 
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the canonical Jordon-Wigner transformation, the Hamiltonian of such a SP chain placed in 
an external magnetic field H can be mapped onto that of interacting (spinless) fermions 
system with band filling tuned by H . An convenient way to deal with ID fermi systems is 
bosonization. After adopting the bosonization technique, the Hamiltonian of such a ID fermi 
system is transformed into a modified sine-Gordon model with a linear spatial modulation 
in the cosine potential. In the boson picture, the effective Hamiltonian is 



Here 4>(x) is a boson field operator, H(x) is its conjugate momentum, U(x) = — *^tjj, and 
they satisfy the commutation relation 



A is a spatial modulation parameter, which represents the effect of fermi surface shifting 
from half filling in the fermion picture. A 

Although the Zeeman energy influences the comparison between the ground state energy 
of boson models with different modulation, it has no effect to the structure of the ground 
state for fixed band filling case. It is neglected in present paper. 

An alternative realization of the Hamiltonian (1) with fixed A refers to the interaction 
fermion system with tunable particle concentration. For example, the charge-density-wave 
degree freedom of the ID extended Hubbard model is described by this model. Away 
from the half filling, or A 7^ 0, people usually assumes that the cosine potential can be 
neglected due to the rapid oscillation. However, according to this statement, the energy 
gap would drop suddenly to zero at A = 0. We believe it is more likely that the cosine 
potential might experience an intermediate lattice phase before entering into the free boson 
field regime. 

In absence of A, the Hamiltonian is a standard quantum sine-Gordon model. Many works 
have been done about this case both in field theory and in condensed matter physics. UHll 
By a variational method Coleman @ first discovered that the energy density of the system 




(1) 




(2) 
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is unbounded below when the coupling constant [3 2 exceeds 8-/T, this corresponds to the 
Kosterlitz-Thouless (K-T) phase transition by the equivalence between the 2D Coulomb gas 
and sine-Gordon model, the ground state is a spatially homogenous commensurate phase. 

In presence of finite A, the spatial modulation parameter favors the IC phase whereas 
the cosine term itself favors the commensurate phase, this makes the model more compli- 
cated. Previously, many results for the Hamiltonian (1) have been obtained in the classical 
limit. 0~@ In the quantum case, this cosine potential is usually considered to shrink to 
zero with any finite A, the crossover between commensurate and IC phases would happen 
abruptly. Schulz has discussed a similar quantum system as Hamiltonian (1) @, he found 
that the domain wall structure appeared in the ground state in the case of the chemical po- 
tential exceeding a threshold value and the energy gap was finite, for large chemical potential 
the presence of the gap was unimportant. In fact, his result implied that there existed a 
soliton lattice phase which had finite density of domain walls between commensurate phase 
with finite gap and IC phase with vanishing gap. 

After taking a shift <p(x) + ^§ — > (j)(x), the Hamiltonian (1) seems as same as that of 
Schulz. But if one pays attention to the boundary condition of boson operator <p(x), he will 
find that two models are different since the boundary conditions of system will be altered 
under this shift, namely, J L d x <f)(x)dx = — > J L d x (p(x)dx = therefore the physical 
properties will be also changed. We prefer the former (zero charge sector), since this integral 
is proporational to the particle number after normal ordering in the fermion picture. To 
our knowledge, with finite A, the complicated configuration between commensurate and IC 
phases has not been explored clearly. 

In present paper, we develop a new method to give a quantitative discussion on the 
soliton lattice phase of the Hamiltonian (1). Our basic idea refers to a variational procedure 
with respect to a spatially inhomogenious Gaussian wave functional which is introduced to 
determine an approximate ground state of the model. Our results show that as long as A 
shifts from zero, the soliton lattice phase appears in the ground state for f3 2 < Air. Only 
in the regime of (3 2 > An the IC phase with vanishing gap can emerge in the ground state 
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for small a. The paper is set up as following. In section II, we introduce a new spatially 
inhomogenious Gaussian wave functional and calculate the energy and the renormalization 
mass of the ground state. The phase diagram and the evolvement of the soliton lattice 
phase with varying spatial modulation parameter are given in section III. Section IV is our 
conclusion. 

II. SPATIALLY INHOMOGENIOUS GAUSSIAN WAVE FUNCTIONAL 

METHOD 

In the Schrodinger field picture, the problem is to solve the functional Schrodinger equa- 
tion 

= (3) 

where \1/ is a wave functional of the variable (f>(x). It is obviously difficult and far beyond 
our present abilities to solve this equation for the Hamiltonian (1). Here we try a variational 
method by choosing some general ansatz for \l>, then we calculate the expectation value of 
the Hamiltonian and minimize the energy with respect to the variational parameters in the 
ansatz. 

If we consider the field operator 4>(x) = 0, the cosine term will oscillate in the space yield- 
ing zero contribution to the ground state. This may not be the most favorable configuration 
of the ground states. One may shift (f>(x) to compensate Ax as much as possible to lower 
the value of the cosine potential term, however, the integration of d x <f> should be zero since 
it must satisfy the boundary condition. Thus (f)(x) should take a kink after some distance, 
moreover, the modulation of (p(x) induces a "distortion energy" \(d x (f)) 2 . Nevertheless, we 
expect an appropriate form of (f)(x) which can lower the total energy at most. 

Taking above into consideration, we approach the ground state of the system by a trial 
Gaussian wave functional with a field shift as 
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V(P, (p d , f) = N f exp ji J P(x)(f)(x)dx - l - J dxdy [((f>(x) - <j> d (x)) f(x, y) (4>(y) - My))}} ■ 

(4) 

Where Nf is the normalization factor, P(x), 4>d{x) and f(x, y) are the variational functions. 
The expectation value of the Hamiltonian of Eq.(l) with respect to the wave functional of 
Eq.(4) is given in as 

E(P, <f) cl , f) = JdxS^ [P(x) 2 + {d x 4> cl {x)f] - ^Zoos ((3<f> cl (x) + Ax) 

1 If d 2 ) 

+-f(x,x)--J dy5{x-y)^-f- 1 (x,y)\ , (5) 



where 



Z = exp|-^r 1 (x,x)| , (6) 
f _1 {x,y) denotes the inverse of f(x,y), i.e., 

J f(x,x')f- 1 (x\y)dx' = 5(x-y) . 
For simplicity we choose f(x,y) of the form 

f{x, y) = ^ J dk^J k 2 + \i 2 cos k(x - y) 

with inverse 

where jj 2 is a variational parameter. 

Up to now, we have introduced three variational parameters P(x), <fi c i(x) and y 2 which 
will be determined in the following calculation. For the ground state, the minimal energy 
configuration is clearly achieved with P(x) = 0. 

Introducing a upper cutoff A in the integral of Eq.(7), Eq.(6) is explicitly evaluated as 



§1 

4tt 



.l + V'i + OiA- 1 ; 
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(8) 



AVl + (M- 1 ) 2 1 



a 



+ - (d x (f)ci(x)) - -^Z cos (f3(f) d (x) + Xx) 



(9) 



and the expectation energy is 

E(<p chf i 2 ) = jdx { _ , 2 ^ XT ^ n p 

After minimizing the energy expectation value with respect to n and (fi c i(x) respectively, 
we get 

/i 2 = a %— J dxcos9(x) , 
d 2 9(x 



dx 2 



= aZ sin 9(x) . 



(10) 



Here, we have set 8(x) = (3<f) c i(x) + Xx for brevity, L is the macroscopic length of the system. 
Eq.(10) is a soliton equation, integrating the above equation with variation of 9(x), we 

get 

— = ±J2aZ (c-cos6>) 



where c is a integration constant and it satisfies 

--K 



2nV^Z 



(11) 



(12) 



Vc+1 V Vc + 1 / A 
according to the boundary condition J L d x (j>dx = 0. Here K is the first class elliptical 

integral function. The renormalization mass equation is simplified as 

2y/2c 



H = aZ 



:K 



2^2 {c+l)E 



(13) 



v/cTT \\c + ij v v " ' Wc + i 

E {^J-^ij is the second class elliptical integral function. The approximate energy density of 
the ground state is 



Ayi+(M-T i 



aZc - 2^ 2 - ^A 2 



(14) 



III. PHASE DIAGRAM AND SOLITON EVOLVEMENT OF THE GROUND 

STATE 



Since Eq.(ll) involves elliptical function, it is difficult to be solved explicitly. Before 
solving it we give some analytic discussion at special conditions for c = 1 or oo. 
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A. Two special cases 



i) When c = 1, we know AA 1 = 0. The non-trivial solution of Eq.(ll) represents one 
soliton, 

cos 26 = ± tanh V aZx . 

For this case Nakano and Fukuyama have given a detailed discussion. '0 What we should 
note is the modification of excitation gap near this case. After expanding the elliptical 
functions to linear terms of AA -1 , we get a expression of the renormalization mass as 

(M- 1 ) 2 = (/ioA- 1 ) 2 - . (15) 

7T 

/XqA -1 is the renormalization mass with AA -1 = 0, it satisfies 



SI 
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(^oA- 1 ) 2 = «A- 2 ( pt= ] . ( I G ) 

',1 + v/l + ^oA- 1 ^ 2 



If f3 2 < 8n, the above equation has a non-trivial solution with nonzero value of fi^A^ 1 . 
Eq.(15) implies that the excitation gap will trail off when AA^ 1 shifts from zero. It should 
be noted that the ground state at finite AA -1 , whatever small, differs from the ground state 
at AA -1 = 0, since the ground state is spatially inhomogenious with soliton lattice structure 
in the first case, whereas it is spatially homogenous commensurate phase in the second 
condition. 

ii) When c — > oo, we get 



— / dxcos9(x) = a ^ , 
LJ K J 2(AA 



dxcos9(x) = n ^ x A _ n2 , (17) 



so the renormalization mass equation is simplized as 

si 

[ } { V2XA-^ l v l + ^l + ( M -i) 2 y' 

In the case of f3 2 < Air, Eq.(18) always has a nonzero solution for the renormalization 
mass fiA" 1 . When (3 2 > 4n, existence of a nonzero solution depends on the competition 



between ahr 2 and AA _1 . In the special case of (3 2 = An, it has only zero solution when 
ahr 2 < 2v / 2AA~ 1 , while it has nonzero solution when ahr 2 > 2y / 2AA~ 1 . For this special 
case the model can be solved exactly by mapping it into a modified Thirring model, which 
predicts similar result that the finite mass appears for sufficiently large ah~ 2 . 

B. The phase diagram of the model 

In order to give an understandable phase diagram, we numerically solve the renormal- 
ization mass equation and soliton equation for a pair of parameters (ah~ 2 , j3 2 ). The families 
of curves with constant renormalization masses for different values of AA" 1 are depicted in 
the parameter plane, see Fig.l. 

With any finite AA _1 , we find that the ground state is always in the soliton lattice phase 
with finite renormalization mass in the range of (3 2 < 4tt. The IC phase with vanishing gap 
only appears in the range of (3 2 > 4%, increasing of AA _1 leads to enlarging of IC phase 
area. It is clear that the transition is a continuous transition on the vertical line (3 2 = An 
for ahr 2 < 2v / 2AA~ 1 since the renormalization mass tends to zero from the left side of 
it. The soliton lattice phase with finite renormalization mass is allowed in the range of 
P 2 > 47T. The phase diagram is shown in Fig.2 where the transition boundary consists of the 
envelope of the family yuA" 1 =constant with different values of AA^ 1 . Under the boundary 
the renormalization mass vanishes, where the Hamiltonian (1) reduces to that of a free 
massless boson field, whereas upon the boundary the renormalization mass is nonzero. The 
first order transition occurs along this boundary in the range of [3 2 > 4ir. 

In the area where the renormalization mass is finite, a soliton lattice phase emerges in 
the real space in the ground state. The local structure of the ground state are shown in 
Fig. 3. It can be seen that the amplitudes of the oscillating part of the local solitons become 
smaller and the periods become shorter with increasing AA _1 . In the soliton lattice phase, 
the real space in the ground state will split into many periodic domains 27r ^~ 1 - > < x < 2?r( ^ +1 ) 
(n integer). In the area where the renormalization mass is zero, the amplitude of solitons 
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reduces to zero, the soliton lattice phase disappears, the ground state enters into a free boson 
regime. 



IV. CONCLUSION 

In this paper we have studied the behavior of the ground state of a sine-Gordon model 
with linear spatial modulation. The main task is to calculate the ground state and low 
excitations for constant modulation parameter. We have applied a variational procedure 
with respect to a Gaussian wave functional to find an approximate ground state. Our result 
shows that for f3 2 > An IC phase grows with increasing A, where the low energy physics can 
be described by a free massless boson field theory, while upon the boundary the excitation 
gap is still finite, and the real space is no more homogeneous, a soliton lattice state appears. 
On the other hand, for f3 2 < Atc, only the soliton lattice phase can exist. With increasing 
spatial modulation parameter, the soliton periodic length decreases, while the correlation 
length (inverse of the gap) increases. The transition type between soliton lattice phase and 
IC phase maybe of the continuous or first order, which relies on the varying of parameter 

Applying above conclusion to the CDW degree freedom of the ID extended Hubbard 

2*t-2V § 



model whose parameter (3 = 8iry 2 ^t+u+5V > we ^ nc ^ ^ ne f°U° w i n g assertion. The ground 
state is always in the soliton lattice phase in the case of U + 13V > 6nt, the umklapp term 
can not be neglected for any filling case. On the other hand when U + 13V < Qnt the 
ground state maybe tuned to a free massless boson state when the band filling shifts away 
half filling sufficiently, only in this case the umklapp term can be directly cancelled. 
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Captions 

FIG.l The families of curves with constant renormalization masses in the parameter 
oA." 2 - f3 2 plane with (a) AA" 1 = 10" 6 ,(b) AA" 1 = 10" 3 and (c) AA" 1 = HT 1 . The 
renormalization mass is chosen as /iiA" 1 = 1, /^A -1 = 10 -1 , /^A" 1 = 10~ 3 , /^A" 1 = 10~ 5 , 
^A" 1 = 10~ 7 , ^A' 1 = 10- 9 and ^A" 1 = 10" 11 . 

FIG. 2 Phase diagram of the model in the parameter aA~ 2 — (3 2 plane with different values 
of AA" 1 . The regions I and II represent soliton lattice phase and IC phase respectively. 

FIG. 3 The local forms of the soliton lattice phase for AA" 1 = 0.004 (a) and AA" 1 = 0.01 
with aA" 2 = 0.2 and (3 2 = 5n. Here, we take^^ as function of coordinate x, since it is 
proportional to the effective cosine potential. 
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